Split Instability of a Vortex in an Attractive Bose-Einstein Condensate 
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An attractive Bose-Einstein condensate with a vortex splits into two pieces via the quadrupole 
dynamical instability, which arises at a weaker strength of interaction than the monopole and the 
dipole instabilities. The split pieces subsequently unite to restore the original vortex or collapse. 
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Quantized vortices in gaseous Bose-Einstein conden- 
sates (BECs) offer a visible hallmark of superfluid- 
ity Q ||, [H , where repulsive interatomic interactions play 
a crucial role in the vortex stabilization and lattice for- 
mation. Attractive BECs, on the other hand, cannot 
hold vortices in any thermodynamically stable state. A 
fundamental issue of the decay of a many-particle quan- 
tum system may be addressed if a vortex is created in 
an attractive BEC. Such a state has become possible ow- 
ing to the development of the Feshbach technique Q] by 
which the strength and the sign of interactions can be 
controlled j§ §, §. 

Suppose that a singly quantized vortex is created in 
a BEC with repulsive interaction and that the interac- 
tion is adiabatically changed from repulsive to attrac- 
tive. According to previous work the vortex state 
remains metastable until a dimensionless strength of in- 
teraction g to be defined later reaches a critical value 
9m(< 0), and when \g\ exceeds \g]^\, the system develops 
a monopole (breathing-mode) instability and collapses. 
In this Letter, we show that this conclusion holds only 
when the system has exact axisymmetry, and that even 
an infinitesimal symmetry-breaking perturbation induces 
the quadrupole dynamical instability that appears for \g\ 
smaller than loj^J. We note that similar dynamical insta- 
bilities ||, [Io| |ll| initiate vortex nucleation observed by 
the ENS group [|2j . A dynamical instability is also shown 
to transfer a vortex from one to the other component of 
a binary BEC system Here, we show that yet an- 

other dynamical instability causes a vortex to split into 
two pieces that revolve around the center of the trap. 
Surprisingly, in some parameter regimes, the pieces sub- 
sequently unite to restore the original vortex, and this 
split-merge cycle repeats. We report below the results of 
our studies on the collapsing dynamics of a vortex in an 
attractive BEC. 

We first investigate the Bogoliubov spectrum of a 
single-vortex state. The single-vortex state is deter- 
mined so as to minimize the Gross-Pitaevskii (GP) en- 
ergy functional within the axisymmetric functional space 
'00 = f(f,z)e 1 ^ with r — (x 2 + y 2 ) 1 / 2 , where we ignore 
the effect of vortex bending [l4|, [ll| . In the following 
analysis, we normalize the length, time, and wave func- 
tions in units of d = (Ti/mu)\_) 1 / 2 , wj 1 , and (N/d^) 1 ^ 2 , 
where uj± is the radial trap frequency, and N, the num- 



ber of BEC atoms. We obtain the Bogoliubov spectrum 
by numerically diagonalizing the Bogoliubov-de Gcnncs 
equations |16) 



(K + 2g\ip a \ 2 ) u n + gip 2 v v 



(la) 



(K + 2g\ip \ 2 ) v n + gip* 2 u n = -E n v n , (lb) 

where K = -V 2 /2 + (r 2 + A 2 z 2 )/2 - ^ with A = u z /uj_, 
and n is the index of the eigenmode. Here, g = AirNa/do 
characterizes the strength of interaction, where a is the 
s-wave scattering length. For a vortex state -00 oc e 1 ^ , 
each angular momentum state u n oc e im ^ is coupled only 
to v n oc e l ( m - 2 )<P i an d we shall refer to m as the angular 
momentum of the excitation. 

We find that there is at least one negative eigenvalue 
in the m = mode for any g < and A even in the 
presence of a rotating drive. The vortex state with at- 
tractive interactions is therefore thermodynamically un- 
stable, and eventually decays into the non- vortex ground 
state by dissipating its energy and angular momentum. 
At sufficiently low temperatures in a high- vacuum cham- 
ber, however, the thermodynamic instability is irrelevant, 
since the energy and angular momentum are conserved. 
In fact, recent experiments Q have demonstrated that 
the vortex state in a stationary trap has a lifetime of 
~ 1 s, which is much longer than the characteristic time 
scales of the dynamics that we shall discuss below. 

When the complex eigenvalues emerge in the Bogoli- 
ubov spectrum, the amplitude of the corresponding mode 
grows exponentially in time. As noise is inevitable in 
experimental situations, such dynamical instabilities are 
more important than the thermodynamic one at low tem- 
perature. 

Figure |l| shows the real and imaginary parts of the 
lowest eigenvalues of the m = — 1 and 3 excitations in 
an isotropic trap. The eigenvalues become complex at 
the critical strength of interaction gg = —15.06, showing 
the onset of the dynamical instability in the quadrupole 
mode. The imaginary part of the complex eigenvalue is 
proportional to ^Sq — 9 as shown in the inset in Fig. [l]. 
The complex eigenvalues emerge also in the dipole modes, 
i.e., m = and 2, for g < g^ = —18.02. The eigenval- 
ues with other m are real for g larger than the critical 
value for the monopole (radial-breathing-mode) instabil- 
ity gfj = -23.7. 
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FIG. 1: For g > gg = —15.06, the solid curve denotes the 
eigenvalue (real) ol the m = — 1 mode, and the dashed one 
that of the negative-norm branch of the m — 3 mode in an 
isotropic trap. These two branches merge at g = gg, below 
which the eigenvalue becomes complex and the dynamical in- 
stability sets in. Inset: The solid line shows the imaginary 
part of the complex eigenvalues for g < gg, which is propor- 



tional to ^J~g\ 
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FIG. 2: The A = u) z /uo± dependence of the critical strengths 
of interaction for the quadrupole mode gg, the dipole mode 
g c £, and the monopole mode g% of the single- vortex state and 
the critical strength of interaction for the non-vortex ground 

State ^nonvortex ■ 



Fi gure |2| shows the A dependence of Qq, g^, 9m i 
and flSnvortex in axi-symmetric traps, where s" onvortcx 
is the critical value for the non-vortex state to collapse 
through the monopole instability. We note that \g M \ 
is always larger than \g^\ and and hence the lat- 

ter instabilities arise before the monopole instability sets 
in. For a trap with A > 0.3, the quadrupole instability 
arises before the dipole one, unlike the quasi-lD toroidal 



trap 0, [[I, [19 • 

We also performed numerical diagonalization for 2D 
systems. Strong confinement in the z direction pro- 
duces the quasi-2D trap, when Huj z is much larger than 
the characteristic energy of the dynamics. The effective 
strength of interactions in the quasi-2D oblate trap is 
given by g 2D = ^X/(2Tr)g 3D Q. In the 2D system, the 
dynamical instability arises in the quadrupole, dipole, 
and monopole modes at £f| Dcr = -7.79, gfj 3 "' = -11.48, 
and <?fp cr = —23.4. The dependencies of the complex 
eigenvalues on g are similar to that in Fig. |l|. 

To understand the dynamical instabilities analytically, 
let us consider the GP action integral in 2D 



dt / drip* 



-I- 1 h - -0 V>- 

dt 2 2 2 m ' 



(2) 



We assume that the state 
^ = Em c m (i)0 m (r), where 4> r , 
form of m (r) = [r' m l / {yji:\m 
im(j>} and d = [1 + g/^ir)] 1 ^ 



evolution is described by 
m (r) is assumed to take the 
l+1 )]exp[ 



Trfl™l +1 )]exp[-rV(2d 2 ) 
minimizes the GP energy 
functional for the m = 1 state. Substituting this ip into 
Eq. (g) and minimizing S with respect to c m yields 
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(3) 



mra\ni'2'i'a?, 



where e m = 



J dr(\V^ m \ 2 + r 2 ^ 



| 2 )/2 and G2\'™\ 



When the BEC exists in the m = 1 
ac, wu uuioui u^) = e~ v * + 0(\c m ^i\ 2 ) with /x = 
' 2 + d 2 + g/(Aitd 2 ). The linear analysis of Eq. (|j) for 
= e v *c m (m/1) yields 



(e m - /x)c m + 2G 

771 1 



G" 



m,2 — m ~# 



(4) 



It follows from this that, for m = —1, the eigenfrc- 
quencies are given by A ± yB, where A = [g/(8w) — 
1]/[1 + ^/(STr)] 1 ^ an d B = 3 + 5 9 /(8tt) + [1 + <?/(27r) - 
5 2 /(32tt 2 )]/[1 + g/(87r)]. We find that B is a monoton- 
ically increasing function for g > — 87r, and B becomes 
negative for g < g" ~ —9.2, which is in reasonable agree- 
ment with (7q Dci = —7.79 stated above. We also find that 
the imaginary part appearing for g < g cr is proportional 
to y/g CI — g, in agreement with the inset of Fig. [j]. 

The Bogoliubov analysis described above is valid only 
if deviations from a stationary state are small. To follow 
further evolution of the wave function, we must solve the 
time-dependent GP equation. Since we are studying the 
growth of small perturbations, high precision is required 
in the numerical integration, and hence we consider the 
GP equation in 2D 
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g 2D W 



(5) 



to ensure sufficiently small discretization in the Crank- 
Nicholson scheme [^l). This situation corresponds to an 
oblate trap with large A. 
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FIG. 3: The density and phase (insets) profiles of the 
time evolution of the vortex state. The initial state is the 
stationary solution of the Gross-Pitaevskii equation plus a 
small symmetry-breaking perturbation. The insets present 
the gray-scale plots of the phase from (2n — l)n to (2n + l)ir 
with an integer n. (a)-(f) Time evolution with <? 2D = —9. 
(g) An expanded image at t — 17, where the interaction is 
switched from g 2 ° = —9 to g cxpan d = 50, and the trap is 
switched off at t — 16. (h) An expanded image at t = 18 with 
ffexpand = 0. The sizes of the images are 7 x 7 in (a)-(f) and 
18 x 18 in (g) and (h) in units of (h/muj^) 1 ^ 2 . The sensitiv- 
ity of imaging in (g) and (h) is 20 times higher than that in 
(a)-(f). 



Figures || (a)-(f) depict the time evolution of the 
density and phase profiles with g 2B = —9, which is 
smaller than the critical value for the quadrupole mode 
<?g Dcr = —7.79 but larger than that for the dipole mode 
gfP cr = —11.48. A small symmetry-breaking perturba- 



tion is added to the initial state to imitate noise in re- 
alistic situations. Due to the quadrupole instability, the 
vortex is first stretched [Fig. || (b)], and then splits into 
two clusters that revolve around the center of the trap 
[Fig. H (d)] with angular velocity ~ 0.73wj_. In the first 
deformation process the m = — 1 and 3 components grow 
exponentially, and their Lyapunov exponents agree with 
the imaginary part of the complex eigenvalues. Interest- 
ingly, the split process is reversible: the two clusters sub- 
sequently unite to restore the ring shape [Fig. || (f)], and 
this split-merge process repeats. We numerically checked 
that no split-merge phenomenon occurs for g 2B > <?q Dci , 
where the system is metastable. The insets in Fig. |3| illus- 
trate the phase plots. At t = 16, there are three topolog- 
ical defects: the central one exists from the outset, and 
the other two enter as the vortex splits, in accordance 



with the fact that the m — 3 component grows upon 
the vortex split. The two side vortices cannot be seen in 
the density plot, and hence they may be called "ghost" 
vortices jllj] that carry very little angular momentum. 

The two clusters in Fig. [| (d) may be regarded as re- 
volving "solitons" whose phases differ by tt. In fact, at 
an energy only slightly below that of Fig. || (d), there 
is a low-lying state in which two solitons revolve with- 
out changing their shapes. It is interesting to note that 
this situation is similar to the soliton-train formation 
observed by the Rice group p2[ , where the modulation 
(dynamical) instability causes a quasi-lD condensate to 
split into solitons when the interaction is changed from 
repulsive to attractive using the Feshbach resonance. 
We note that similar instabilities split an optical vortex 
propagating in a nonlinear medium into spiraling soli- 
tons |2^, . This similarity between attractive BECs 
and optical solitons [ p6[ implies that other nonlinear phe- 
nomena, such as pattern formation, which has been pre- 
dicted in attractive BECs p7| , may also be realized in 
optical systems. 

When the system becomes too small to be observed by 
the in situ imaging method due to the attractive interac- 
tion, the condensate must be expanded before imaging. 
Figure || (g) shows the expanded image at t = 17, where 
the interaction is switched from g 2D = — 9 to g exp and = 50 
and the trapping potential is switched off at t = 16 [Fig. || 
(d)]. The image shows the interference fringes due to the 
overlap of the atomic clouds emanating from the two clus- 
ters. The wavelength of the interference pattern is pro- 
portional to the expansion time. Figure |^ (h) shows the 
expanded image at t = 18 with ^expand = 0. Comparing 
Figs. [| (g) and (h), we find that the stronger repulsive in- 
teraction produces more fringes and bends them around 
the center. 

When |# 2D | exceeds |gf) Dcr | = 11.48, a dipole instabil- 
ity arises in addition to the quadrupole one. The dipole 
instability causes atoms to transfer from one cluster to 
the other, thereby inducing the collapse. Figures ^ (a)- 
(c) show the collapse process with g 2D = —11.5. After 
the split-merge process repeats a few times, the balance 
between the two clusters is broken due to the dipole in- 
stability. As a consequence, the cluster labeled A grows 
[Fig. U (a)], then B grows [Fig. |] (b)] like a seesaw, and 
eventually B absorbs most atoms and collapses [Fig. |] 
(c)], where the original topological defect begins to spiral 
out as indicated by the white arrow in the inset. With a 
stronger attractive interaction g 2B = — 12, both clusters 
collapse immediately after the vortex split as shown in 
Figs. H (d)-(f). In this collapse process, we found the ex- 
change of a vortex-antivortex pair (see the insets) . This 
phenomenon is also seen in the split-merge process with 
g 2D = —11.5, while it is not seen at weaker attractive 
interactions, say, at g 2D = — 9. 

Figures [| (g)-(i) show the collapse where the interac- 



tion is switched from g = to g 
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FIG. 4: The density and phase (insets) profiles of the time 
evolution after the vortex split with g 2D = —11.5 in (a)-(c) 
and with g 2D = -12 in (d)-(f). In (g)-(i), g 2 ° is switched 
from to —24 at t — 0. The right inset in (e) illustrates the 
topological defects and branch cuts in the left one. Here we 
note that the end part of the central branch cut is separated 
to yield a vortex-antivortex pair 3-4 [inset of (c)]. Vortex 5 
and antivortex 4 are subsequently combined [inset of (f)]. The 
sizes of the images are 5 x 5 in (a)-(f), 8 x 8 in (g), 4 x 4 in 
(h), and 2 x 2 in (i) in units of (Ti/majx) 1 ' 72 . The sensitivity 
of imaging of the density plots is 1/4 in (a)-(f), 1 in (g), 1/5 
in (h), and 1/20 in (i) in units of those in Figs. N (a)-(f). 



t = 0. The vortex first shrinks due to the monopolc in- 
stability [(g) — ► (h)], then splits into two clusters due to 
the quadrupole instability [(h) — > (i)], and both clusters 
collapse. When the interaction is switched to an even 
greater attractive one, a shell structure is formed [ p7[ , 
which splits into several parts due to multipole instabil- 
ities, and each fragment collapses and explodes pro- 
ducing very complicated collapsing dynamics. 

In conclusion, we have studied the dynamical instabil- 
ities of a quantized vortex in an attractive BEC. The dy- 
namical quadrupole instability spontaneously breaks the 
axisymmetry and splits a vortex into clusters that re- 
volve around the center of the trap, which then unite to 
restore the vortex or eventually collapse. The dynamical 
instabilities presented here play a much larger role than 
the thermodynamic one at low temperature, and serve as 



a dominant mechanism for the collapsing dynamics of a 
rotating condensate: vortices collapse via the dynamical 
instabilities around the topological defects. 
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